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SUMMARY 


The response of a thin, rigid viscoplastic plate subjected to a spatially 
axisymmetric Gaussian step-pressure impulse loading was studied analytically. 

A Gaussian pressure distribution in excess of the collapse load is applied to 
the plate, is held constant for a length of time, and then is suddenly removed. 
The plate deforms with monotonically increasing deflections until the dynamic 
energy is completely dissipated in plastic work. The simply supported plate of 
uniform thickness obeys the von Mises yield criterion and a generalized consti- 
tutive equation for rigid viscoplastic materials. For the small deflection 
bending response of the plate, the governing system of equations is essentially 
nonlinear. Transverse shear stress is neglected in the yield condition and 
rotary inertia in the equations of dynamic equilibrium. A proportional loading 
technique, known to give excellent approximations of the exact solution for the 
uniform load case, was used to linearize the problem and obtain the analytical 
solutions in the form of eigenvalue expansions. The linearized governing equa- 
tion required the knowledge of the collapse load of the corresponding static 
problem. 

The effects of load concentration, an order of magnitude change in the 
viscosity of the plate material, and load duration were examined while holding 
the total impulse constant. In general, as the load became more concentrated, 
both the peak central velocity and the time for plate motion to cease increased. 
For the less viscous plate, these increases of velocity and time were more pro- 
nounced. The final plate profile became more conical as the load concentration 
increased but did not approach the purely conical shape predicted for the ideal 
point impulse by the rigid perfectly plastic analysis with the Tresca yield cri- 
teria. Profiles of the less viscous plate were influenced more by the load con- 
centration. The shorter the time of load application, the greater the initial 
kinetic energy, the plate response, and the final deformation. 


INTRODUCTION 

There have been a number of plastic circular plate analyses (refs. 1, 2, 
and 3) which permit a time variation of the load; however, there have been few 
papers which consider a radial variation other than linear (refs. 3 and 4). 
Florence's solutions (ref. 5) for uniform pressure acting over a central region 
permit one to vary a single parameter to obtain all intermediate states between 
the point load and uniform distribution. The only other general spatial distri- 
bution of load which has received significant analytical attention is the Gaus- 
sian distribution. By varying a single parameter, this general distribution can 
span the extremes from the point load to the uniformly distributed load. This 
versatility was recognized by Sneddon (ref. 6) who approximated the dynamic load- 
ing of a projectile on a thin, infinite elastic plate by a Gaussian distribution 
of pressure. Madden (ref. 7), in his study of shielding of space vehicle struc- 
tures against meteoroid penetration, related the meteoroid-shield debris loading 



of the main vehicle wall to a Gaussicin initial velocity distribution. The first 
study of this loading on a plastic plate was made by Thomson (ref. 8). He 
obtained the solution for a rigid, perfectly plastic plate obeying the Tresca 
yield condition subjected to an ideal impulse of Gaussian distribution. The 
corresponding solution for the rigid viscoplastic plate obeying the von Mises 
yield condition is presented in reference 9. Weidman (ref. 2), in considering 
the response of simply supported circular plastic plates to distributed time- 
varying loadings, presented an example of a radial Gaussian distribution of 
pressure with an exponential decay. The plate material in reference 2 was also 
rigid, perfectly plastic obeying the Tresca yield condition. An approximate 
solution for a uniformly distributed load acting over a circular region of 
radius a, assuming a Gaussian distribution of velocity at all time, was 
obtained by Lepik (ref. 10). 

The present paper presents the analytical results for the small deflection 
bending response of a simply supported circular plate of rigid viscoplastic 
material. The plate is subjected to a spatially axisymmetric Gaussian distri- 
bution of pressure having a rectangular pulse in time. (The corresponding solu- 
tion for a Gaussian ideal impulse loading is presented in ref. 9.) The effects 
of load concentration, of an order of magnitude change in the viscosity of the 
plate material, and of load duration are examined while holding the total impulse 
constant. Approximate expressions are developed for the time at which plate 
motion ceases, for the final shape of the plate, and for the final central 
displacement . 

Part of the information presented in this report was included in a thesis 
entitled "Gaussian Impulsive Loading of Rigid Viscoplastic Plates" submitted by 
Robert J. Hayduk in partial fulfillment of the requirements for the Degree of 
Doctor of Philosophy in Engineering Mechanics, Virginia Polytechnic Institute 
and State University, Blacksburg, Virginia, August 1972. 


SYMBOLS 

Aji series coefficients, equation (A6) 

Ci,C 2 constants defined by equations (A5) 


F 

f 


nondimensional load amplitude, 
yield function 


\/3PoR^ 

4Mo 


h plate half-thickness 


Jp(x),Ip(x) Bessel function of first kind of real and imaginary arguments, 
respectively 


J 2 second invariant of deviatoric stress tensor 

kp,K(j, radial and circumferential curvature rates 


2 



k 

Mo 

m 

n 

P» 

TiMo 


Po 

Po 

Q 

q 

R 


r 





t 


tf 

U(P) 

u(p,t) 


V 


w 

z 

a 

3 

Y 


yield stress in simple shear 
yield moment of plate, 

radial and circumferential bending-moment resultants 

radial and circumferential bending-moment resultants at initial yield 

nondimensional radial bending-moment resultant, M^/Mq 

nondimensional circumferential bending-moment resultant, ^^/MQ 

P^r2/i _ e-g> 


total nondimensional load-carrying capacity. 


Mo \ 


3 


pressure amplitude at plate center 

nondimensional pressure amplitude at plate center, PqR^/Mo 
shear stress resultant 

nondimensional shear stress resultant, RQ/Mq 

plate radius 

radial coordinate 

deviator stress tensor 

deviator stress tensor at initial yield 

time 

time for motion to cease 

steady component of velocity 

dynamic component of velocity 

nondimensional plate velocity, i”. 

aMo 3t 

transverse deflection of plate 
transverse coordinate 
plate geometry and material constant, 

nondimensional Gaussian shape parameter 
material constant , y°/2k 


\/3liR^Y 

2hM„ 


, sec 


3 



$(f) 

<KXn,3) 

4'(An,P,3) 

V2,V^ 


material constant 
strain rate tensor 

eigenvalues determined from equation (A7) 
mass density per unit area of plate 
nondimensional radial coordinate, r/R 
stress tensor 

yield stress in simple tension 
load duration, sec 
function defined by equations (3) 
circumferential coordinate 
function defined by equation (A13) 

function defined by equation (A12) 
harmonic and biharmonic operators in cylindrical coordinates 


LINEARIZATION OF THE GENERALIZED CONSTITUTIVE EQUATIONS 
FOR RIGID VISCOPLASTIC MATERIALS 

A generalized constitutive equation for rigid viscoplastic materials is pre 
sented in this section. Material elasticity is neglected in order to simplify 
the analysis as is frequently done in theoretical investigations of dynamic plas 
tic response of structures. Such an assumption is generally believed to be 
valid when the dynamic energy is considerably larger than the maximum energy 
which could be absorbed in a wholly elastic manner and the duration of loading 
is short compared with the fundamental period of vibration. 

Perzyna (ref. 11) developed a constitutive equation for rate-sensitive plas 
tic materials by generalizing the constitutive relationships used by previous 
researchers (Hohenemser and Prager in ref. 12; and Prager in ref. 13). The gen- 
eralized constitutive equation proposed by Perzyna is 

= Y° $(f) (1) 

aaij 

where is the strain rate tensor and 





k 


1 


( 2 ) 



The second invariant of the deviatoric stress tensor is 





where 


Sij = c^ij - ^ 5ij<^kk 


The function in the generalized constitutive equation (eq. (1)) is defined as 


$(f) = 0 
$(f) i 0 



(3) 


Also, k is the yield stress in simple shear and y® denotes a physical con- 
stant of the material. 


Perzyna (ref. 14) has shown that the generalized constitutive equation for 
viscoplastic materials reduces to the constitutive equations of an incompressi- 
ble, perfectly plastic material first considered by von Mises and to the flow 
law of perfect plasticity theory. As in the theory of perfectly plastic solids, 
convexity of the subsequent dynamic loading surfaces and orthogonality of the 
inelastic strain-rate vector to the yield surface follow from Drucker's postu- 
lates defining a stable inelastic material with inclusion of time-dependent 
terms (ref. 11). 


A method of linearizing boundary-value problems in the theory of visco- 
plastic solids is described by Wiezbicki in reference 15. In this method, as 
shown graphically in figure 1 , the concept of proportional loading is used to 
relate the state of stress in the subsequent stress state to the state 

of stress on the initial yield surface f = 0. For proportional loading 

the direction cosine tensor of the state of stress in deviatoric space is inde- 
pendent of time. The proportional loading that satisfies equation (2) is given 
by 




'ij 


"ij 


(4) 


This equation is a reasonable approximation for axisymmetrically loaded, simply 
supported circular plates because the plate center amd boundary are automati- 
cally proportionally loaded; that is, the bending moments must always be equal 
at the plate center amd the radial bending moment must always be zero at the 
plate boundary. 

When equation (4) is utilized, the generalized constitutive equation 
(eq. (1)) becomes 

Eij = 1 9(f) 2ij (5) 
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s • • 

vrtiere the viscosity constant y equals y°/2k and f = - 1 . The ratio 

Sij/Sij may be expressed in terms of any one of the deviatoric stress compo- 
nents. The simplified constitutive equation (eq. (5)) still is nonlinear in 
stresses. However, for this analysis, the linear form 


$(f) = f 


( 6 ) 


is chosen to produce full linearization of the constitutive equation (eq. (5)) 
as follows : 


'ij 


= X( j - ^ ) 


ij' 


(7) 


It should be noted that equation (7) is really a flow relation for a given struc- 
ture rather than a constitutive equation describing a given material (ref. 16). 

For the problem of a uniformly loaded, simply supported circular plate with 
$(f) B f, Wierzbicki (ref. 15) has shown that the approximate solution, obtained 
with the proportional loading hypothesis, agrees very well with a numerical 
finite-difference solution without the proportional loading approximation. The 
solution of the linearized problem also agrees well with experimental data on 
impulsively loaded plates obtained by Florence (ref. 17). 


GOVERNING EQUATIONS AND BOUNDARY AND INITIAL CONDITIONS 

A Gaussian transverse load, described by 

p(r) = (8) 

is suddenly applied to the entire surface of a rigid viscoplastic plate of 
radius R and thickness 2h, where 3 is a nondimensional load shape parame- 
ter. The boundary of the plate at r = R is simply supported. After being 
applied steadily during the time interval 0 1 t 5 x, the load is suddenly 
removed. The geometry of the plate and loading are shown in figure 2. For 
3=0, equation (8) describes a uniform load; as 3 approaches * and Pg 
approaches ", equation (8) describes a point load at the plate center. 

In the derivation that follows, the equations apply to the loading phase 
0 "I t 1 X . To obtain the equations for the postload phase t > x , the 
applied load terms must be omitted. The internal forces and moments acting on 
a typical plate element are shown in figure 3- If rotary inertia is neglected 
but transverse inertia is taken into account, the equations of motion are 

^(rQ) + rpoe-8r-^/R^ 

3r 

l-(rMp) - Hj, = rQ 
3r 
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With the usual assumptions of small deflection plate theory - that is, 
plane sections remain plane, strains vary linearly through the thickness, and 
strains are linearly related to the curvature - the curvature-rate — moment rela- 
tions, derived from the linearized constitutive equation (eq. (7)), are 




- «p) - <2^4, - Hp)] 


( 10 ) 


The moments Mj, and satisfy, for any value of r, the equation of the 

initial yield surface, 


R2 - + m| = M§ 


( 11 ) 


where Mq = Ooh^ is the yield moment of the plate material and a© = \/3k is 
the yield stress in simple tension. 


For small deflections of the plate, the curvature rates Kj 
related to the deflection rate w by 

a^w^ 


and are 


K. 




K - ^ 

4) “ r 3r 


( 12 ) 


Equations (9), (10), and (12) form a linear parabolic system of partial differ- 
ential equations with six unknown functions - M^, Q, w, K^, and ^ - 

plus the unknown static moment distributions and^ M^. 

By eliminating all unknowns except w, the system of governing equations 
can be reduced to the single, fourth-order equation: 


^Mph 

3\l3y 


V^w - Ppe“^^^/^^ + ji = 1 ^ 


3t r 3i 


l-(rM^) - 


(13) 


where 


= 

32 + 1 9 

32 + 1 9 


3 r 2 r 3r 

_9r2 r 3r 


The right-hand side of equation (13) represents the internal force distribution 
in the plate at the initiation of collapse for the static case. 

Let P(J denote the amplitude of the static load-carrying capacity of the 

plate, then the right-hand side of equation (13) can be replaced by -p^e'^r^/R^ 
and the governing equation becomes 
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/ 


V^w + y M = (pq - p^)e-3r2/R2 

3t 

This method of solution, proposed by Wierzblcki (ref. 15), has the impor- 
tant property of replacing the unknown static moment distributions and 5i(j,, 

whose explicit formulas are not known for the von Mises yield condition, by the 

load distribution at initiation of static collapse, that is, Thus, 

the need for explicit formulas has been reduced to finding the value of a con- 
stant p^ corresponding to a particular value of the shape parameter B. The 
determination of the load-carrying capacity of a circular plate under a Gaussian 
distribution of pressure is presented in reference 18. For completeness the 
variation of Pq and total nondimensional load-carrying capacity, 

P» _ - e~^ \ 

ttMq Mo \ e / 


^Mph 

3\f3y 


obtained by integrating the pressure distribution over the plate, are shown as 
functions of 3 in figure 4. 


Define the nondimensional quantities as 


Mr 

m = _£ 
Mn 


n 


M. 


- 

M„ 


P = £ 
R 


(15) 


PR2 3w 
aMp 


F = 


_ \/3PoR2 


, . \/3pAR2 


4m; 


F' = 


"5m; 


and let + 1 then the final form of the governing equation (eq. (14)) 

3p2 P 9p 

is 


V^v + 3 a ly = 2\/3(F - F')e-3P^ 
2 3t 


( 16 ) 


where 


a = N/3yyR^ 


2hM^ 


The boundary conditions of the simply supported plate are 

m = n and q n 0 (p = 0) 

v=m=0 (p=1) 


(17) 


By using equations (10), (12), and (9), equations (17), in terms of rate 
of deflection, become 


lim/3£v _ J 3v\ - 0 
p->-oUp2 P 3pj ■ 


(18a) 
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lim/ 33y ^ 2 3^v _ 1_ _ q 

P'^0\3p3 P 3p2 p2 3p/ 



2 + I = 0 

9p2 P p=1 ' 


(18c) 


v(1,t) =0 (I8d) 

For the Gaussian step-pressure loading, the plate is initially flat and at rest; 
that is , 

w(p,0) = 0 v(p,0) = 0 (19) 


RESULTS AND DISCUSSION 


The solution to the governing equation (eq. (16)) with associated boundary 
and initial conditions (eqs. (18) and (19)) is presented in the appendix. Effects 
of load distribution and plate viscosity on plate response are examined while 
holding the load duration t and the total impulse I^ = ttR^tPo3“Vi - e“6) 
constant. The constants used in the graphical presentations were chosen to 
be consistent with results of reference 15; that is, x □ 0.1 msec and 






B 10. Values of a of 1 x 10“3 and 1 


sec are 


used to show effects of plate material viscosity. Effects of load distribution 
are given by the parameter g. The required variation of Pq with g in 
order to maintain constant total impulse is given by the parameter 




B 

1 - e-2 


The values of 


determined from equation (20) and used in the computations are 
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The load becomes more concentrated at the center of the plate as 0 is increased 
and the amplitude increases almost linearly as 0 becomes large. The graphical 
results were obtained by programing the modal solution as given in the appendix 
(eqs. (A14), (A15), (A17), and (Al8)). 

In figure 5, the time histories of the plate central velocity for 0 rang- 
ing from 10“3 to 10 000 are shown collectively for the two values of the parame- 
ter a. The effect of varying the shape of the load from a uniform one to essen- 
tially a point load is seen immediately. The peak velocity, which occurs at 
the instant of load removal, is increased tremendously - a factor of approxi- 
mately 6 for a = 1 X 10“3 sec and 10 for a = 1 x 10“^ sec. The time for the 
central velocity to become zero also increased with 0, this effect being more 
pronounced for a = 1 x 10~2 sec. 


The results of figure 5 indicate smaller values of nondimensional velocity 
for a = 1 X 10~2 sec. However, the nondimensional velocity v(0,t) plotted 
is the actual velocity 3w/3t divided by the constant aMg/pR^. Thus, although 
the nondimensional value of velocity decreases with an increase in a, the 
actual velocity increases. For example, for 0 = 10 000 and t n t, the curves 
of figures 5(a) and 5(b) show that 


which corresponds to 


= 0.291 


|H(0,t)' 

3t 


a=1x10 


-2 


i"(0,T) 

3t 


= 2.91 


-‘a=1x10-3 


Thus, the less viscous plate (the one with the higher value of a and, conse- 
quently, the higher value of y) attains higher velocities, has larger displace- 
ments , and acquires more energy from the applied load which requires longer times 
to dissipate than the more viscous plate. 


The profiles of the plate are shown in figure 6 for both values of a and 
various values of 0 . As the load becomes more concentrated the profile becomes 
more conical, as one would expect. Figure 6(b) indicates that the profiles of 
the less viscous plate are influenced more by the shape parameter 0 than are 
those for a = 1 x 10”3 sec. 


An approximation of the deflection of the plate is obtained from equa- 
tion (Al8) by retaining only the first term of the series and by using the 
approximation 
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CO 

Cl + C2p2 + 1. e-3p^ ■*■(-+ P^] (-1)”(3p^)” =, -1^ J_ 4)(Xi ,p,B) 

2B U / ^ (2n)n! 3 ^ ’ 


which is obtained from equation (A14) with v(p,0) = 0. The result for 
T < t < tf is 


Mfi- w(p,t) 

Qt^Mo 


= ^ 
4Ba 


Cl + C2p2 + 1_ e-Bp2 + / 1 + p2\ ^ (-1)n(Bp2)" 

2B U ) ^ (2n)n! 


X (fF't - Ft - Ja 


(F 


_ p,)g-(2XT/3a)t 


_ pg“(2X 1 /3ot) ( t-T ) ^ p, 


( 21 ) 


An approximate expression for the time for motion to cease can be obtained by- 
setting the derivative of the approximate displacement expression to zero, that 

is, ^1 = 0, is 

tf 


at 


tf □ 2®- In 


2x!( 


1 + 


F (p(2X-i/3a)T 


( 22 ) 


Equation (22), plotted in figure 7 for 0 1 B 1 100 and the two values of a, 
is an implicit function of B since F and F' vary with B* The effect of 
B , as shown in figure 7 , diminishes after an initial rapid rise of tf with 
increasing B • The symbols represent computed times using the complete equation 
for the velocity (eq, (A17)). Equation (22) is a very good approximation for 
a = 1 X 10“3 sec; however, except for small values of B» the approximation is 
poor for a = 1 X 10”2 sec. 


The substitution of equation (22) for tf into equation (21) provides an 
approximate expression for the final plate displacements 


w(p,tf) = \/l_ 
a2M, ^ 


Cl + C2p2 + :L e-3p2 + /I + p2\ (-1)»(Bp^)" 

^ ^ 2B U / ^ (2n)n! 


X (¥' 3o_ In 

2X^ 


1 + |l(e^2Xi/3a)T _ ^ 
TT« ' 


- Ft 


and an approximate expression for the final center displacement 


w(0,tf) = ^/Ci + LUF' ^ In 

a2M, 2 bA 2x4 


, , f^(e< 24 / 3 a)T . ,) 

F« 


- Ft 


(23) 


(24) 
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Equation (24) is plotted as a function of 3 for the two values of ot in fig- 
ure 8. The approximations are in excellent agreement with the points computed 
from the exact equations for both a = 1 x 10”3 sec and 1 x 10 "^ sec, although 
the tf approximations for the larger value of ot were poor for large values 
of 3 as shown in figure 7. Again, the nondimensional central displacements 
are shown smaller for a = 1 x 10“2 sec whereas the actual displacements are 
larger than for a = 1 x 10”3 sec. 

Although not shown, profiles obtained from the approximation (eq. (23)) 
were compared with profiles obtained from the exact equation (eq. (Al8)). For 
a = 1 X 10~3 sec, the differences between the approximate and exact profiles 
were negligibly small for the entire range of 3 considered, 10”3 to 10 000. 

However, for the less viscous plates, a = 1 x 10 “^ sec, the differences were 

not negligible and the approximation (eq. (23)) should therefore be restricted 
accordingly. 

The effect of load duration on the final center displacement is shown in 
figure 9. Total impulse was held constant as before. For these results the 
approximate expression (eq. (24)) was used with a n 1 x 10~3 sec and 3 = 10. 
The point T = 0 represents the ideal impulse and was calculated by using the 
corresponding approximate center displacement equation of reference 9. As seen 
in the figure, the center displacement monotonically decreases as the pulse dura- 
tion lengthens. The maximum center displacement occurs for the ideal impulse. 

The final center displacement becomes zero at T = 3-32 x 10“^ sec, the impulse 

duration at which the applied load becomes equal to the collapse load of the 

plate. For impulses beyond that duration, the plate remains rigid. 


CONCLUDING REMARKS 

A thin, simply supported, rigid, viscoplastic plate subjected to a Gaussian 
step-pressure impulse loading has been analyzed within the realm of small deflec- 
tion bending theory. The plate material obeys the von Mises yield criteria and 
a generalized constitutive equation. These considerations lead, essentially, to 
nonlinear equations governing the dynamic response of the thin plate. A propor- 
tional loading hypothesis, known to give excellent approximations of the exact 
solution for the uniform load case, was used to linearize the problem and obtain 
analytical solutions in the form of eigenvalue expansions. The linearized gov- 
erning equation on the velocity of the plate required the knowledge of the col- 
lapse load of the corresponding static problem, that is, the collapse load for 
the specific load distribution parameter 3. 

The effects of load concentration, an order of magnitude change in the 
viscosity of the plate material, and load duration were examined while holding 
the total impulse constant. In general, as the load became more concentrated, 
both the peak central velocity and the time for plate motion to cease increased. 
For the less viscous plate material, these increases of velocity and time for 
plate motion to cease were more pronounced. The final plate profile became more 
conical as the load concentration increased but did not approach the purely coni- 
cal shape predicted by the rigid, perfectly plastic analysis with the Tresca 
yield condition for an ideal point impulse. As the viscosity of the plate 
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decreases, the shape parameter has more effect on the final deformed plate 
profiles. 

Approximate expressions were developed for the time at which plate motion 
ceases, the final shape of the plate, and the final central displacement. Com- 
parisons with the exact series solution indicated that all approximations were 
excellent for a plate property constant a of 1 x 10“3 sec throughout the range 
of shape parameter 0. For a = 1 x 10 “^ gee, the approximation for the final 
central deflection was good for the entire range of 3, but the other approxima- 
tions were limited in usefulness. 

Comparisons between the results of the step-pressure loading and the ideal 
impulse loading indicate that the viscoplastic plate subjected to a certain 
total impulse achieves a larger peak velocity and, consequently, kinetic energy 
if the impulse is delivered instantaneously - that is, as an ideal impulse - 
rather than over a finite length of time. Consequently, the viscoplastic plate 
subjected to the ideal impulse requires a longer time after load removal to dis- 
sipate its higher energy and experiences larger permanent deflections than the 
plate subjected to the step-pressure loading of the same total impulse. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 8, 1978 


APPENDIX 


SOLUTION OF GOVERNING EQUATION BY EIGENVALUE EXPANSION 

The governing equation (eq. (16)) can be solved by means of an eigen expan- 
sion method since the right-hand side of the equation is not a function of time. 
Equation (16) is repeated here for convenience: 

V^v + 3 a ^ = 2\j3(F - F')e-®P^ 

2 at 


Substitution of 

v(p,t) = u(p,t) + U(p) 
into equation (16) results in 

V^u(p,t) + 3 a gu (P>t) + v%(p) = 2J3(F - F')e~^P^ 

2 at 


(A1) 


which separates into 

V^u + 3 a lli = 0 

2 at 


(A2) 


and 


V^U = 2\/3(F - F')e-^P^ 


(A3) 


Equation (A3) is the same as equation (16) except for the absence of the inertia 
term. Thus, U(p) is an equilibrium solution of equation (16) with the same 
boundary conditions (eqs. (18)), as follows: 


U(p) 


= -V1(F _ 

43 


C2P‘ 


23 


.-3p2 




z 

m=1 


(-i)°^(3p^ 

(2m)m! 


where 


Cl 


C2 


L - 1 - L e-^ - 1 (-1)”g°^ 

63 6 33 3 “(2m)m! 

m= I 


> 

OO 

-1. + I - 1_ e-2 - 

63 6 63 (2m)m! 

m= 1 . 


(A4) 


(A5) 


A general 
all prescribed 


solution due to Wierzbicki (ref. 15) satisfying equation (A2) and 
boundary conditions can be written in the form 
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u(p,t) =2An[lo(Xn) Jo(XnP) " «Jo(Xn) ^ (A6) 

n=1 

where Jq(x) and Iq(x) denote the Bessel functions of the first kind of real 
and imaginary arguments. The solution (eq. (A6)) identically satisfies the 
boundary conditions (eqs. (I8a), (l8b), and (l8d)). The eigenvalues are 

roots of the following transcendental equation stemming from the boundary con- 
dition (eq. (I8c)) of zero bending moment at the plate edge: 

lo(Xn) Jl(Xn) + Il(Xn) Jo(Xn) - 4X„ lo(Xn) Jo(Xn) = 0 (A7) 

The only remaining unknowns in the solution are the series coefficients 
A^. These coefficients are evaluated from the initial condition (eqs. (19)), 
that is, 

v(p,0) = u(p,0) + U(p) =0 

Thus, 

u(p,0) = -U(p) (A8) 

and substituting equation (A6) for u(p,0) gives the following equation: 

00 

An[lo(Xn) Jo(XnP) - Jo(Xn) lo(XnP)] = -U(p) (A9) 

n=1 

The coefficients A^ can be determined from equation (A9) by virtue of the 
orthogonality of the system [lo(Xn) Jo(XnP) - JQ^Xn) Ifl^XnP)] on the interval 
[]0,l] where p is used as a weighting function. Therefore, coefficients A^j 
can be determined as 


An = - 


r pU(p) I^Io(Xn) Jo(XnP) - JQ^Xn) ^O^XnP)] dp 

•'o 

f p[lo(Xn 
Jc\ 


(A10) 


) Jo(XnP) - Jo(Xn) lo(XnP)!^ dp 


where U(p) is defined by equation (A4). The resulting coefficients are 

An[lo(Xn) Jo(Xnp) - Jo(Xn) lo(XnP)] = - F' )1_ ilj(Xn,P ,8) (All) 

3 


with the functions \l)(Xn»P,8) defined by the relation 

-IM iKXn,p,3) = <|)(Xn,8)|'lo(Xn) Jo(XnP) - Jo(Xn) lo(XnP)] 


(A12) 



t 


APPENDIX 


where 

and 


N 




(A13) 


N = 1 Cl lo(Xn) Jl(^n) - ^ Ci JqCXh) + C 2 Io(Xn)rj- 


^n 




■-!» Jl(Xn) +-^ Jo(Xn) 




- C 2 Jo(Xn) 


1 Il(Xn) -2 lo(Xn) +i_ Ii(Xn) 
Xn X2 X3 


+ L. lo^^n) f xe"^^^ Jo(XnX) dx 
2 g Jq 


- i' 


xe“3x2 Iq(x^x) dx 


°o^ 

0 ’ m=1 


(-1)°>(3x2) 

(2m)m! 


m 


JoCXj^x) dx - Jq(X 


"> I is 


\ 00 

(-1)“(3x2)®‘ 

m= 1 

L (2m)ml 


lo(XnX) dx 


Phase 1 : 0 < t < t 

When equations (A4) and (A6) are summed and equation (All) is used, the 
velocity during loading, phase 1 (0 < t < t), becomes 


v(p,t) = - il(F - F') 
43 




(-1)»^(3p2)'^ 

(2n)n! 



1 _ 


\|)(Xn,p ,3 )e 


-(2X|*/3a)t 


(A14) 


The displacement of the plate is determined by integrating equation (A14) 
with respect to time. Taking the initial condition of zero displacement 
(eqs. ( 19 )) into account, the displacement becomes 
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w(p ,t) 

a^Mo 


_^(F - F') 
43a 


Cl + 


Cop^ + 1— e”^P^ 
23 



(-1)m(3p2)°» 

(2m)m! 


t 



4^ >P ^ 


(e-(2X^/3a)t 



(A15) 


Phase 2: x < t < tf 


After the load is removed from the plate at time 
the governing differential equation becomes 


T , phase 2 (x < t < tf ) , 


V^*v + 3 a ^ = -2\/3F'e-3P^ 

2 at 


(A16) 


Utilizing the same eigenvalue expansion techniques as before with continuity of 
velocity from phase 1 to phase 2 produces the plate velocity during phase 2 
motion as follows: 


v(p,t) = l/3 F’ 

43 


Cl + C2p2 


I_ e"2P^ 
23 



(-1)n(3p2)» 

(2n)n! 


^ ^ \p(Xn>P>8) 


(F - F')e‘ 


-(2\^/3a)x _ p g-(2Xf*/3a)(t-x) 


(A17) 


Integration of equation (A17) with respect to time with continuity of dis- 
placement to determine the arbitrary function of integration produces the dis- 
placement of the plate during phase 2 motion as follows; 


pr2 

a^Mo 


w(p,t) = ^/LCF't - Fx) 



+ C?p2 + J_ e”^P^ 
23 


(1 + p^\ Y' (-i)°^(3p^)°* 

^ (2m)m! 


+ 


2 \/ 3 X' i_ 4»(Xn.P,3) 

^rfx9 



p, yg-(2X^/3a)t 


pg-(2X^i/3a)(t-x) ^ p, 
(A18) 


Equations (A14), (A15), (A17), and (Al8) represent the complete solution 
for the velocity and displacement of the plate. In the limit as 3 approaches 
0, the Gaussian distribution becomes the uniform distribution of load and this 
solution reduces to the solution presented by Wierzbicki (ref. 15). 
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Figure 1.- Representation of proportional loading in deviator ic stress space 



Z 


Figure 2.- Simply supported circular plate with Gaussian distribution 
of pressure applied during time interval T. 
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Figure 3.- Element of circular plate with applied forces and moments. 



Figure 4.- Nondimensional collapse load amplitude and total load-carrying 
capacity as functions of load distribution parameter 3. 







Figure 7-- Comparison of approximate expression (eq. (22)) for time for 
motion to cease tf and points determined from complete equations for 
step-pressure loading. 



0 

Figure 8.- Comparison of approximate expression (eq. (24)) for final central 
deflection and points determined from complete equations for step-pressure 
loading. 


24 



T, msec 

Figure 9-- Variation of final center displacement with load duration 
impulse held constant; a = 1 x 10 “3 sec; 3 =10. 
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